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Abstract 

In this work, the effect of introducing next nearest neighbor (NNN) hopping to the 2D materials was 
studied using the graphene 2D honeycomb two sublattice as an example. It is found that introducing NNN 
to the 2D materials removes the symmetry around the Fermi level and shifts it, at a small value of NNN 
hopping. This effect increases with increasing NNN hopping. If the NNN hopping becomes competitive 
with nearest neighbor (NN) hopping, the dispersion relations of the 2D materials changes completely from 
NN hopping dispersion relations. The results show that the 2D material sensitivity for NNN hopping effect 
is much larger in the 2D honeycomb lattice than 2D square lattice. This is due to the fact that the number 
of NNN sites is equal to six, which is the double of NN sites in the 2D honeycomb lattice. Meanwhile, the 
number of NNN sites is equal to four which is equal to NN sites in 2D square lattice. We therefore conclude 
that by changing the ratio between NNN and NN sites in the 2D lattice one can tune the sensitivity for NNN 
hopping effects. 



* malimed62@uwo.ca 



1 



I. INTRODUCTION 

It is shown in [[T]-31| that the physical properties of 2D materials like graphene and magnetic 
stripes are mainly attributed to both their lattice structure and range of interactions between its 
sites. With fixing the range of interactions to include only the nearest neighbor hopping, a com- 
parison between the obtained results of the 2D square lattice, zigzag edged, and armchair edged 2D 
honeycomb lattice show that the 2D lattice structure and its edge configuration play very important 
rule in its dispersion relations and consequently its possible applications. 

Despite this importance of lattice structure in the properties of 2D materials, experiments and 
theories [l5]-0 show that increasing the range of the interaction to include the next nearest neigh- 
bor (NNN) in the graphene 2D honeycomb lattice changing its dispersion relations by removing 
dispersion symmetry around the Fermi level with shifting Fermi level value, and changing their 
behavior around the impurities in the lattice [|9||. Also, including the next nearest neighbor hopping 
in 2D square lattice changing its dispersion relations [lOJ . 

It is interesting to study the effect of introducing the next nearest neighbor in the structure of the 
hopping E matrix, and consequently on the obtained dispersion relations and the localized edge 
states of the 2D materials which has not been previously examined. 

In this work, the graphene Hamiltonian [6] which includes the next nearest neighbor interaction 
term will be used to study its eff'ect on the dispersion relations, edge states, and impurities states 
in the graphene nanoribbons with zigzag and armchair edge. The obtained results should also be 
applicable to the magnetic case, since the next nearest neighbor interaction term can be added to 
Heisenberg Hamiltonian [fm[T2ll . 

n. THEORETICAL MODEL 

The system initially under study is a 2D graphene nanoribbon in the xy-plane. The crystal- 
lographic description of graphene honeycomb lattice is given in Q. The nanoribbon is of finite 
width in the y direction with atomic rows (labeled as r = 1, • ■ ■ , A'^) and it is infinite in x direction 
(-00 ^ 00). The total Hamiltonian of the system is given by 




(1) 
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FIG. 1 . Armchair (left) and zigzag (right) graphene 2D honeycomb nanoribbons in xy-plane, where black 
(gray) dots are the sublattice A(B) with a line of impurities (white dots) in the middle of the sheet, with 
average spin alignment in z direction. The stripes are finite in y direction with A'^ rows {n = 1 , • • • ,N) and 
they are infinite in the x direction. Figure taken from llT3ll . 

where the first term ?o,;(~ 2.8eV) is the nearest- neighbor hopping energy and in graphene it is the 
hopping between different sublattices A and B. Also O.leV) is the next nearest-neighbor 
hopping energy which here in graphene is the hopping in the same sublattice ^W\. The sum- 
mations over i and j run over all the sites where i and j are belong to different sublattice A(B) 
for the nearest neighbors hopping term, and they are belonging to the same sublattice for the next 
nearest-neighbor hopping energy. Where the nearest neighbor hopping to.j has a constant "bulk" 
value t when either i and j are in the interior of the nanoribbon, and another constant value when 
i and j are both at the edge of the nanoribbon (i.e., in row n = I or n = N). Similarly, for the next 
nearest-neighbor hopping energy t^.., we assume that it has a constant value f when the site / is 
inside the nanoribbon, and it is equal to t'^ for sites at the edge of the nanoribbon. 

Since the nanoribbon extends to +00 in the x direction, we may introduce a ID Fourier trans- 
form to wavevector along the x direction for the fermions operators aj (a,) and b^j {bj) as follows: 



Here A^o is the (macroscopically large) number of spin sites in any row, Qx is a wavevector in the 
first Brillouin zone of the reciprocal lattice and both r, and is the position vectors of any carbon 




(2) 
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sites / and j. The new operators obey the following commutation relations: 



(3) 



Also, we define the hopping sum: 



V 



■>-.,) 



-(q.v-(r,-r,) 



(4) 



The sum for the hopping terms ?o/i,y is taken to be over all v nearest neighbors and over all v' 
next nearest-neighbor in the lattice which depends on the edge configuration as zigzag or armchair 
for the stripe. For the armchair configuration, the hopping sum for nearest neighbors gives the 
following factors T„n'iqx) 



Tnn'iqx) — t 

and for the zigzag configuration, it gives 



n',n±l 



(5) 



'^nn'i.^x) ^ 



2 COS 



V3 



-q^a 



^n',n±l + ^n',n+l 



The hopping sum for next nearest neighbors gives the following factors t^„/(<?a) 

Tnn'iqx) = t' [5n',n±2 + 2 CO?,{qxai ll)6n' ,n±\\ 

for the armchair configuration, and 

Tnn'iqx) = 2t' [cos( ^f3qxa)6n',n + cos( ^qxal2)6n',n±2[ 



(6) 



(7) 



(8) 



for the zigzag configuration case, where the + sign, in all the above factors, depends on the sub- 
lattice since the atom line alternates from A and B. 

Substituting Equations ([2]) and (|4]) in Equation ([1]), and rewriting the summation over nearest 
and next nearest neighbors sites, we get the following form of the operator term H: 



H = -Y_i ^'^^^^ i^n^^"' + ^nbn') + T(qx)anbl, + T{-q;)albn' 



(9) 



nn' 



The first terms count the elementary excitations on each sublattice, while the second describes the 
coupling between the sublattices. 

In order to diagonalize H and obtain the dispersion relations for graphene nanoribbons, we 
may consider the time evolution of the creation and the annihilation operators a] {ai) and b^- (bj), as 



calculated in the Heisenberg picture in quantum mechanics where the time dependent is transferred 
from the system wavefunction to the operators. In this case, the equations of motion (using the 
units with /i = 1) for the annihilation operators ai{bj) are as follows [fT4l - [T8l : 

ddfi 



dt 



i[H, an] 



(10) 



and 



db„ 

--i=i[H,b„] 
dt 



nn' 



(11) 



where the commutation relation in Equation Q was used, as well as the operator identity {AB, C] = 
A[5,C] + [A,C]5. 

The electronic dispersion relations of the graphene (i.e., energy or frequency versus wavevec- 
tor) can now be obtained by solving the above operator equations of motion. The electronic energy 
can be expressed in terms of the frequency using the relation E = fioj, and assuming that electronic 
energy modes behave like e,x^{-ico{qx)t], on substituting this time dependent in Equations and , 
we get the following sets of coupled equations: 



oj(q^)an = ^ Tnn'(qx)a„' + T„„^(-qx)bn' 

nn' 

Oj(qx)bn = ^ Tnn'(qx)an' + T'„n'iqx)bn' 
nn' 

The above equations can be written in matrix form as following 

coiqx) 

where the solution of this matrix equation is given by the condition 



(12) 
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T'iqx) T(q,) 


a,n 


(14) 
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T*(qx) T'iq,) _ 


bn 





det 



-{(x){q^)lN - T'{q^)) T{q^) 

riq,) -ioj(qx)lN - T'iq,)) 



(15) 



Where T^q^) and T'{qx) are nearest and next nearest exchange matrices respectably, which are 
depend on the orientation of the ribbon, and co{qx) are the energies of the modes. The matrix Tiq^) 



TABLE I. Nearest neighbor hopping matrix elements for the graphene as 2D honeycomb lattice 

Parameter Zigzag Armchair 
a te-^'i-" 

y t te'i^all 

TABLE n. Next nearest neighbor hopping matrix elements for the graphene as 2D honeycomb lattice 



Parameter Zigzag Parameter Armchair 

e It'cosi^qxO) 6 t' 

^ 2f co?,{^qxal2) ri 2?'cosfea3/2) 



is given by 



or ;S 

/3 a r 

J a J3 

/3 a 



The matrix T'(qx) for zigzag ribbon is given by 

^60^00 
6 ^0 

^ 6 ^ 

^ 6 
0^06 



and the matrix T'{qx) for armchair ribbon is given by 

^ T] 9 

r] T] e 
6 T] r] 

e T] J] 
e T] 



(16) 



(17) 



(18) 



The parameters 7, e, ^, and 77 depend on the stripe edge geometry and are given in Tables |I] 
andim 



A. Neglecting the next nearest neighbor hopping as special case 



The next nearest neighbor hopping t' can be neglected compared to nearest neighbor hopping 
t, in this case the T'{qx) is equal to zero matrix and Equation ([14]) become as following 



det 







(19) 



-{oj{qx)lN) T{q;) 
T*{q^) -(a){q_^)lN) 

which is the result obtained before for graphene ribbons using the tight binding model with ne- 
glecting NNN hopping [031 . It is also very similar to the case of magnetic stripes in (21 , which 
do not have NNN exchange, the only difference between the magnetic case and TBM graphene 
without NNN hopping is the effect of a, i.e. insite energy. 



III. NUMERICAL CALCULATIONS 

The dispersion relations for the above graphene nanoribbons are obtained numerically as the 
eigenvalues |fT9ll20l for the matrix Equation This is very similar to graphene with only (NN) 
[|2l, and therefore the same numerical calculations method used there will be used here to get its 
solutions. 



IV. RESULTS 

To compare our results for NN and NNN interaction with the tight-binding Hamiltonian, with 
only NN interactions, we choose our stripes sizes, scaling our result to be dimensional less quan- 
tities, and choose physical parameters matched that ones used in reference lfT3]| for graphene. 

Figures |2] shows the effect of next nearest neighbor interaction in the dispersion relations, edge 
states, and impurities states in the graphene zigzag nanoribbons, as expected all Figure show the 
removing dispersion symmetry around the Fermi level with shifting the Fermi level value toward 
valance band. Figure[2](a) and (b) show the effect for the next nearest neighbor value of t' = 0.036?, 
which correspond to f « 0.1 eV and ? « 2.8 eV given in references [l6l|7]|> The changing in the 
dispersion symmetry around the Fermi level with shifting the Fermi level value toward valance 
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FIG. 2. The effect of next nearest neighbor interaction in the dispersion relations, edge states, and impurities 
states in the graphene zigzag nanoribbons. Right side stripe width N = 20 (a) t' = 0.036f (b) t' = OAt (c) 
t' - 0. 1? and with impurities hne at row number 1 1 with Jj - 0. Left side stripe width N -21(d) t' - 0.036? 
(e) t' - O.lt ({)f- OAt and with impurities Une at row number 11 with Jj - 0. 
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FIG. 3. The effect of next nearest neighbor interaction in the dispersion relations, band gap, and impurities 
states in the graphene armchair nanoribbons. Right side t' - OAt for stripe width (a) N - 20 (h) N - 21 
(c) = 22. Left side f - OAt and with impurities line at row number 1 1 with Jj - for stripe width (a) 
A^ = 20(b)A^ = 21 (c)A^ = 22. 



band is small compared by the obtained results in ||2l for the same zigzag nanoribbons without 
NNN interaction, as the value of NNN increases to f' = OAt as shown in Figure [2] (c) and (d) 
the density of states increase in the conduction band (high energy levels) and decreasing in the 
valance band (low energy levels) removing the symmetry around the Fermi level and shifting 
it. The Figures show that including NNN effecting the flatness of the edge localized states of 
zigzag graphene nanoribbons but not effecting its position in Fermi Level, as the NNN increase 
the flatness decreases which reflect the introducing of depends for hopping in edge sites, which 
more clear for the extended localized edge state in zigzag with width N = 21. 

Figures [2] (e) and (f) show the effect of NNN on the impurities states in the zigzag nanoribbons. 
It is clear that the position of energy state of impurities line not affected by including NNN, which 
is a result of not participating for the impurities in NNN hopping in this calculation. But introduc- 
ing the NNN hopping in the lattice shifting the Fermi level and changing the density of the states 
around the impurities level. As, NNN increase the impurities level move to more density of states 
region, this explain the appearance of impurities level as a moving peak in the density of states for 
the graphene, with increasing surrounding density of states as NNN increasing IfTOl . 

Figures |3] show the effect of next nearest neighbor interaction in the dispersion relations, and 
impurities states in the graphene armchair nanoribbons. The behavior in the armchair case with 
including NNN is very similar to the zigzag case given above for removing the symmetry around 
the Fermi level, shifting it, and its effect in the impurities level relative position to Fermi level. 
There is no any effect on the absence of edge states in armchair nanoribbons and on their band 
shape at qxa/n = with including NNN in the model. 

V. DISCUSSION AND CONCLUSIONS 

In this work, the effect of introducing NNN hopping to the 2D materials was studied using the 
graphene 2D honeycomb two sublattice as example. Including the NNN in the model add NNN 
hopping matrix T'{qx), which depending on the momentum q^ in the direction of nanoribbons 
symmetry, to the diagonal sub matrices uIn in the E matrix. This shows that NNN hopping matrix 
T'{qx) is real and is describing the hopping with translation motion in the sublattice sites in both 
zigzag and armchair stripes. 

When t' is equal to very small percentage of t, the probability is very small for the particles 
to hopping in the same sublattice by NNN hopping and consequently the net number of particles 
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FIG. 4. The effect of high next nearest neighbor interaction t' = 0.5t in the dispersion relations on zigzag 
left side and armchair right side lattice with width A'^ = 20. 

in NNN hopping is very small which result a small effect in the dispersion relations of only NN 
hopping. As the percentage increases, the probability for NNN hopping increases, and conse- 
quently the net number of particles in NNN hopping increases. This results in an increasing effect 
in the dispersion relations of only NN hopping. The main effect of NNN hopping in small range 
is changing the density of states for dominated NN hopping dispersion relations, which can be 
explained as following: since the probability for the particles to hopping in the same sublattice 
by NNN hopping is increasing with increasing its energy and consequently the net number of 
particles from every mode that able to do NNN hopping is proportional to the mode energy. The 
highest energy mode has the highest percentage number of particles that participating in NNN 
hopping, this percentage of particles decreases with decreasing the energy of the mode, most of 
this NNN particles will be trapped in low energy modes. The overall effect is the available parti- 
cles densities is lowest in high energy modes and highest in low energy modes. This means that 
the available momentum spaces in high energy modes is increased for particles in NN hopping 
due to NNN effect consequently the density of states is increased in high energy modes, while the 
available momentum spaces in low energy modes is decreased for particles in NN hopping due to 
NNN effect consequently the density of states is decreased in low energy modes. This removing 
the symmetry around the Fermi level and shifting it, this effect increases with increasing NNN 
hopping. If the NNN hopping become competitive with NN hopping the dispersion will changing 
complectly as seen in Figure |4j 

The above results for NNN hopping is applied to 2D square lattices as shown in Figure |5] by 
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FIG. 5. The effect of next nearest neighbor interaction in the dispersion relations on magnetic 2D square 
lattice with width = 8. Right side t' = OAt left side t' = 0.5t. 

adding the term (S J' /2)(2cos{qxa)) in upper and lower off diagonal of its E matrix. The result 
show the same behavior for the density of states. 

The comparison between the results of f = O.lt between 2D honeycomb lattice and 2D square 
lattice show that the sensitivity for NNN hopping effect is much larger in the 2D honeycomb lattice 
than 2D square lattice, this due to the fact that the number of NNN sites is equal to six which is the 
double of NN sites in the 2D honeycomb lattice, while the number of NNN sites is equal to four 
which is equal to NN sites in 2D square lattice. Therefore by changing the ratio between NNN 
and NN sites in the 2D lattice one can tune the sensitivity for NNN hopping effects. 
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